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Abstract. We present a set of evolutionary computations for binaries with a 
B-type primary at birth. Some liberal computations including loss of mass and 
angular momentum during binary evolution are added to an extensive grid of 
conservative calculations. Our computations are compared statistically to the 
observed distributions of orbital periods and mass ratios of Algols. Conservative 
• Roche Lobe Over Flow (RLOF) reproduces the observed distribution of orbital 

' periods decently but fails to explain the observed mass ratios in the range G 

[0.4-1]. In order to obtain a better fit the binaries have to lose a significant 
amount of matter, without transferring too much angular momentum. 
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■ 1. Introduction 
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lEggletonl ((2000) introduced the denomination "liberal" to make a distinction 



between binary evolution with mass and angular momentum loss and the con- 
^ ■ servative ca se where no mass an d consequently no angular momentum leave 

the system. iRefsdal et alJ (|l974l ) showed that the binary AS Eri is the result 
of liberal binary evolution; the amou nts of mass a nd angular momentum lost 
by the system are however uncertain. ISarnal (|l99.t l showed that only 60 % of 
the mass lost by the loser of (3 Per was captured by the gainer and that 30 % 
of the initial angular momentum was lost during RLOF. Hence it is clear that 
the liberal evolutionary scenario is important for binary evolution calculations. 
However the amount of mass and angular momentum that has to be removed 
$_i ' from interacting systems is far from obvious. 

With the Brussels si multaneous evolut i on co de (a description is given in 
"The Brightest Binaries", IVanb everen et alJ JUIi)) we calculated a representa- 
tive grid of conservative evolu tion of binari es with a B type primary at birth. 
Application of the criterion of iPetersI ( 200 lT ) allows then to determine for each 



of the evolutionary sequences the beginning and ending of various Algol-stages. 
This criterion states that in the semi-detached system: 



The less massive star fills its Roche lobe (RL) 

The most massive star does not fill its RL and is still on the main sequence 
The less massive star is the coolest, the faintest and the largest 



The grid allows to determine an expected distribution of orbital periods and 
mass ratios o f Algols. These can b e compared to the well established observed 

distributions (Budd ing et al. (200l)). 

Considering the examples mentioned above ( Refsdal et alJ ( 1974h . ISarnal 



1993)) it may be expected that the match between the observations and the 
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conservative results are far from satisfactory. Therefore we add liberal calcula- 
tions to our conservative library. The detailed evolutionary tracks can be found 
in http : //www . vub . ac . be/astrof ys/ 



2. Details of the conservative calculation 

Our grid contains binaries with sufficiently small initial periods to lead to Case 
A RLOF (RLOF A): i.e. during H core burning of the donor (initially primary 
star that becomes the less massive after Algol-ignition). The track across the 
HRD is traced for every system i n our grid. A binary lives its era of " Arg olism" 
(|De Loore and Van Rensbergen! (|2004f )) when it obeys the criterion of iPetersI 
(2001). Every binary shows its Algol A (Algol during H core burning) aspect 
for some time during RLOF A. The drastical change of the mass ratio and the 
orbital period during this process can be followed in detail for every system in our 
grid. It happens frequently that RLOF A is succeeded by Case B RLOF (RLOF 
B): i.e. during H shell burning of the donor. Systems that have sufficiently large 
initial orbital periods so that RLOF A does not occur will also show a short 
living Algol B (Algol during H shell burning) appearance. T hese systems have 
been considered in a previous paper (|Van Rensbergen! (j2003)). Figure ^ shows 
a typical example of a case AB. The {7+A.2)Mq binary with an intial period of 
2.5 d starts as Algol A with q=l and P=2.06 d. It ends this Algol A stage after 
some 20 million years with q=0.27 and P=6.87 d. From these inital values the 
system remains an Algol B for some 1.5 million years during H shell burning of 
the donor. The system eventually evolves into a long periodic (P>50 d) Algol 
B with q ~ 0.1. The grid contains some 240 more evolutionary tracks. 
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Figure 1. Conservative evolution of a (7+4-2)Mq binary with an intial pe- 
riod of 2.5 d. The system remains Algol A for ps 20 million years before it is 
Algol B for « 1.5 million years. 
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3. Conservative simululation 

From our grid of conservative calculations we made a simulation of the distribu- 
tion of orbital periods and mass ratios of Algols using a Monte Carlo algorithm. 
As starting conditions we selected: 



• The primaries IMF of lSaloeterl (|1955T ): C(M) + M~ 2 - 35 

• an initial distrubution of orbital periods from lPopova et all (1982): 

n(P)n-i 

• an initial mass ratio distribution as derived bv lVan Rensbergenl ( 200 if ) from 
the non-evolved systems of the catalogue of Spectroscopic Binaries shown on 
http : / / sb9 . astro . ulb . ac . be/ 



The distribution of the initial mass ratio q 



obeys relation (1): 



^(q) (1 + 9) a ; a = 3.37 for early B & 1.47 for late B primaries (1) 

Figure El shows the observed distribution of orbital periods of Algols. Here 
we find more A than B-cases, because the Algol A-phase goes on for a fraction 
of the nuclear time scale, whereas the Algol B-phase lasts only for a fraction of 
the much shorter Kelvin-Helmholtz time scale. RLOF A produces Algols that 
follow the observed distribution rather well. Since the periods of the Algol B 
cases peak towards the long periods, a contribution of a few % of B cases to the 
Algol population will mimic the observed period distribution well. 
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Figure 2. Observed distribution of the orbital periods of Algols compared to 
conservative evolution. Cases B produce an excessive amount of long periods. 
Cases A follow the observed distribution better. There is no doubt that there 
are far more Algol A than Algol B cases. 



Figure shows the observed mass ratio distribution of Algols. The mass 
ratio q is now defined as the observers do: q = ^ donor . 
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The mass ratios of the Algol B cases peak towards the smallest mass ratios, 
whereas RLOF A produces a majority of Algols in the q-bin [0.2-0.4]. As a con- 
sequence an admixture of cases A and B will never reproduce well the observed 
mass ratio distribution in the q-bins [0.4-1] . We may conclude that liberal binary 
evolution is needed to describe the mass ratio distribution of Algols. 
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Figure 3. Observed distribution of mass ratios of Algols compared to con- 
servative evolution. Cases B produce more than 80% Algols with the smallest 
mass ratio q £ [0-0.2]. Cases A produce most of their Algols with q £ [0.2-0.4]- 
The fact that almost 80% Algols are observed with q £ [0-4-1] excludes con- 
servative evolution as the major channel through which Algols can be formed. 



4. Confining the liberal model 

Mass loss is defined by a parameter (3 describing the fraction of the mass lost 
by the loser that is accreted by the gainer: 

M gainer = -(3 M dmor With < < 1 (2) 

Conservative evolution is described with (3 = 1, whereas the liberal case 
uses values of (3 < 1 which are not known beforehand. Conservative evolution 
implies that no angular momentum can be lost by the system, whereas the 
amount of loss of angular momentum in the liberal case is also a free parameter 
if no physics restricts the assumptions. It is clear that realistic hydro dynamical 
calculations should learn us the appropriate choice of the amounts of mass and 
angular momentum which are lost by a binary at any moment of the RLOF 
process. In the mean time we have performed a number of liberal evolutionary 
calculations and compared the results with the observations. 

Our calculations reveal that the time dependent parameter f3(t) should suf- 
ficiently often differ drastically from 1. 

For a given value of (3, the amount of angular momentum lost by the system 
is defined by the position of the site where matter is recoiled from the system. 
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ISoberman et alJ (11997ft have argued that matter can be trapped in a Keplerian 
ring after transit across the second Lagrangian point L<i- The radius of the ring 
is T) times the semi major axis of the orbit. This Keplerian ring passes outside 
L2 which is located at 77 ~ 1.25 . This yields a minim um value of 77 ~ 1.25. 
Hydrodynamical calculations of iLubow and "Shul l|l975ft locate the ring at 77 m 
3. This yields a maximum value of 77 ~ 3. A Keplerian ring located at 77 ~ 2.25 
takes away as much angular momentum as the co-rotating point L2 would do. 
A value of 77 ~ 2.25 is thus a fair value in the interval [1.25-3] that can be used 
to calculate the change of the orbital period of a binary as a consequence of loss 
of mass and angular momentum through a ring which rotates wi t h Ke plerian 
velocity around the center of mass of the system ((Soberman et al.l (11997^ . 

Our calculations reveal that if the system loses angular momentum across 
points located near 77 ~ 2.25 the orbital periods shrink drastically. Most binaries 
become mergers before they show Algolism. The obtained distribution of orbital 
periods of Algols is completely shifted towards the shortest periods. Since this 
conclusion conflicts with observations, the mass that leaves the system carries 
rather the angular momentum of points located near 77 ~ than near 77 ~ 2.25. 



AMI 



□ Beta =1 / 7Mo-primary 
■ Beta=0 .5/ 7Mo— primary 



[< 0] [0-0.2] [0 . 2-0 . 4] [0 . 4-0 . 6 ] [0 . 6-0 . 8 ] [0.8-1] [0-1.2] [ 1 . 2- 1 . 4 ] [ 1 . 4 - 1 . 6 ] [> 1.6] 

Log P 



Figure 4. Liberal evolution with a lot of mass loss ((3=0.5) and a little loss 
of angular momentum leaves the distribution of orbital periods of Algols almost 
unaltered. This figure shows the result for Algols issued from RLOF A and a 
7 Mq primary at birth. 



For all these (statistically sound but physically not well understood) reasons 
we performed a number of calculations with constant (5 = 0.5 and angular mo- 
mentum lost at the edge of the gainer. Until now only the representative cases 
with a 7Mq primary at birth and initial periods leading to RLOF A have been 
performed. Figure 0] shows that the assumption of mass loss through a point 
near 77 = does not alter the conservative distribution of the orbital periods of 
Algols. A conclusion that meets the observations since the conservative orbital 
period distribution matched the observations fairly well. Figure [5] shows that 
our liberal assumption (/3=0.5, r? ~ 0) deviates radically from the conservative 
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case. The q-bins [0.4-1] which were not populated in the conservative scenario 
are now populated properly as required by the observations (see figure 0). 
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Figure 5. Liberal evolution with a lot of mass loss (f3 — 0.5J and a little loss 
of angular momentum changes the mass ratio distribution of Algols drastically. 
The observed mass ratios in the q-interval [ 0.4-1 ] which were not produced 
by conservative evolution are created by this liberal model. This plot shows the 
result for Algols issued from RLOF A and a 7Mq primary at birth. 
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